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Abstract 

The uniaxial structure of Si and Ge (001) facets leads to nontrivial topo- 
logical properties of steps and hence to interesting equilibrium phase tran- 
sitions. The disordered flat phase and the preroughening transition can be 
stabilized without the need for step-step interactions. A model describing 
this is studied numerically by transfer matrix type finite-size-scaling of inter- 
face free energies. Its phase diagram contains a flat, rough, and disordered 
flat phase, separated by roughening and preroughening transition lines. Our 
estimate for the location of the multicritical point where the preroughening 
line merges with the roughening line, predicts that Si and Ge (001) undergo 
preroughening induced simultaneous deconstruction transitions. 
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I. INTRODUCTION 



The structure of the (001) facets of Si and Ge is very interesting from the point of view 
of equilibrium phase transitions. These surfaces have an uniaxial reconstruction [|l|, where 
the uniaxial direction switches by 90° at alternating surface heights. Due to this, the mono- 
atomic and bi-atomic steps have nontrivial topological properties. This atomic structure and 
the crossover from mono-atomic steps in non-vicinal surfaces to bi-atomic steps in vicinal 
surfaces have been studied extensively [0-0]. Roughening type phase transitions in these 
surfaces close to the melting temperature, are an another interesting topic P,[7J. One of 
us suggested earlier f7j that this unusual topology leads to disordered flat (DOF) phases 
and preroughening (PR) transitions without a need for step-step interactions. In this paper 
we present a detailed numerical transfer matrix finite-size-scaling (FSS) study of the model 
introduced in Ref. 0. 

Consider a surface like Si (001), but one which does not reconstruct. Such a surface is still 
uniaxial and it still switches by 90° at alternating surface heights. At finite temperature 
T, thermodynamically-excited steps appear. They separate domains of flat regions. The 
uniaxial structure leads to two distinct types of mono-atomic steps, labeled by Sa and Sb- 
The subscripts denote whether the uniaxial direction in the upper terrace near the step is 
parallel (A) or normal (B) to the step edge. Considering the fact that the uniaxial direction 
switches by 90° at alternating surface heights, one finds, as shown in Fig. 1, that the steps 
have the following topological properties @: (i) If two neighboring parallel steps are of the 
same type, one must be an up-step and the other a down-step, (ii) If a step turns over 90° 
it must change its type, from Sa to Sb and vice versa. Bi-atomic steps exist as well ||, but 
they are probably free-energetically unfavorable close to the roughening temperature @. 

These topological properties imply that terrace excitations have an ellipsoid shape, and 
that the long axes of nested terraces are parallel (perpendicular) if the height change is 
up-down or down-up (up-up or down-down) across the nested terraces. This creates an 
entropic penalty against forming hills and valleys. In other words, it opens the possibility 
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to stabilize a DOF phase without the need for step-step interactions 0. Thus far step-step 
interactions were believed to be crucial for the existence of DOF phases. In this surface 
topology, however, the DOF phase originates directly from the uniaxial structure of the 
surface. 

The restricted solid-on-solid (RSOS) model on a square lattice with the Hamiltonian 
H = E { K (H* + 1,2/)- h(x, y)f + A sin ^(h(x, y) + h(x + 1, y)) } 

+ J2 [K{h{x, y + l)- h{x, y)f - A sin ~(h(x, y) + h{x, y + 1)) } , (1) 

was introduced in Ref. [[/J to describe the thermodynamic properties of such steps in more 
detail. h(r) is an integer- valued height variable at each site r = (x,y). Height differences 
between the nearest neighbor sites are restricted to and ±1. This means that only mono- 
atomic steps are allowed. Bi-atomic steps can be included in a later stage if experimental 
evidence shows they remain important close to roughening temperatures. The model of 
Eq. (|1|) contains two parameters. The A terms distinguish between Sa and S^-type steps: 
Ea = K — A and Eb = K + A are the step energies. Without loss of generality, the uniaxial 
direction is taken to run vertically (horizontally) at even (odd) heights. 

The model Hamiltonian contains two limiting cases. The conventional RSOS model at 
A = displays a Kosterlitz-Thouless (KT) type roughening transition between the flat and 
rough phases ||. On the other hand, in the limit where Ea = and T = 0, Sa steps cost no 
energy while Sb steps are frozen out. In a typical configuration the surface contains a set 
of randomly-placed parallel Sa steps in the form of straight lines. However, the topological 
rule (i) requires that they are alternating up and down steps. This is a typical morphology 
of surfaces in the DOF phase 0. The DOF phase is an intermediate phase between the flat 
and rough phases, where the steps are disordered positionally but have long-range up-down- 
up-down order. It was argued in Ref. 0] that this DOF type structure is stable at finite 
temperatures, in terms of a Fermionic type perturbation theory. 

In this paper, we investigate the phase diagram quantitatively through a detailed transfer 
matrix FSS study. It is important to confirm the existence of the DOF phase numerically. 



The analysis in Ref. f?j was mostly qualitative. The other purpose of this work is to obtain 
a good estimate for the critical value of the ratio r = Eb/Ea, below which the DOF 
phase disappears, see Fig. 2. In real surfaces this ratio takes specific values. For example, 
observations of step fluctuations in STM and LEED experiments yield for Si (001) that 
r ~ 3 (Ea = 325-fT and Eb = 1045/T) [TIJ. By comparing this ratio with the critical value 
r c , one can decide which path the Si (001) surface follows. 

Our model, Eq. ([]]), does not incorporate the 2x1 type dimerized surface reconstruction 
of Si and Ge (001). Therefore it does not describe the competition between surface recon- 
struction and surface roughening in those surfaces. This issue was addressed in Ref. JF|. 
The preroughening line in Fig. 2 is most likely replaced by a PR induced simultaneous de- 
construction transition and the roughening line segment at r < r c by a roughening induced 
simultaneous deconstruction transition. A proper quantitative description of this requires 
at least a RSOS model coupled to an Ising model. We did not study such a model, since the 
number of degrees of freedom becomes too large to obtain meaningful transfer matrix FSS 
results. The precise location of r c in Fig. 2 is the result of a delicate entropy balancing act 
of nested terraces associated with the peculiar 90° switching in the uniaxial direction. Our 
value of r c should be meaningful for Si and Ge (001) if the coupling with the Ising degrees 
of freedom does not change the value of r c by too much, which is a reasonable assumption. 

In Sec. m, we introduce various kinds of interface free energies. They decompose into 
the free energies of Sa and Sb type steps, and show distinct FSS behaviors in the flat, 
DOF, and rough phases. We obtain the phase diagram Fig. 2, by evaluating these interface 
free energies using the transfer matrix method. The numerical results and a summary are 



presented in Sec. II I 



II. INTERFACES AND THE TRANSFER MATRIX FORMALISM 

Consider the model given by Eq. (|1]) on a finite N x M lattice with periodic boundary 
conditions (PBC's), h(x + N,y) = h(x,y) and h(x,y + M) = h(x,y). The ordered flat phase 
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is commensurate with PBC's. Other (gauge invariant) boundary conditions (BC's) create 
frustrations, and thus impose steps in the surface. The interface free energy rj is defined as 
the excess free energy per unit length for each type of BC compared to that of PBC's. Their 
FSS behaviors are different in the various phases. We obtain the structure of the phase 
transitions by studying suitable ones. 

The Hamiltonian in Eq. (fj) is invariant under the global transformations 

h(r) -> h(r) + 2n (2) 
h{r) -> -h{r) + 1 (mod 2) , (3) 

for all integers n. So it is natural to consider the following boundary conditions: step type 
BC's with h(x + N,y) = h(x, y) + 2 and h(x, y + M) = h(x, y); and anti-periodic type BC's 
with h(x + N,y) = —h(x, y) + 1 (mod 2) and h(x, y + M) = h(x, y). We will refer to them 
as HI and H2 respectively. Similarly, VI and V2 refer to the same BC's but with the roles 
of the M and iV interchanged. The interface free energies are defined by 

Va = -±-\n^- , (a = Hl, H2) 
Vp = -hn^- , (0 = VI, V2) . 

iV ^PBC 

with Z a the partition function satisfying the boundary condition a, and all energies and free 
energies measured in units of /c^T. 

Figure 3 shows the topological frustrations induced by these BC's. HI and VI require 
at least two parallel steps; one is an Sa type step and the other an Sb type step (see Figs. 3 
(a) and (c)). Therefore, r/m and 7/vi decompose into t]a + t]b', with tja and T]b the Sa and 
Sb step free energies. On the other hand, H2 and V2 can be satisfied by configurations with 
only one Sa type step (see Figs. 3 (b) and (d)). Therefore, 7] m and r/v2 are equal into t\a- 

These interface free energies must behave in a specific way in each type of phase. The 
step free energy tjb is finite in the flat phase and also in the DOF phase, but vanishes in 
the rough phase. The step free energy tja, is finite in the flat phase, but vanishes in both 
the DOF and rough phase. Therefore, in the flat phase, all four rj a (a=Hl, H2, VI, V2) are 



finite. In the DOF phase, rjm and rjyi remain finite, but rjm and 7/V2 vanish (exponentially 
with system size). In the rough phase, all four rj a vanish (as a powerlaw in the infinite M 
and/or N limit). 

The rough phase is a critical phase where its critical fluctuations are described by the 
Gaussian model. The height- difference correlation function diverges logarithmically with 
distance: 



((h(r) - MO) 2 ) ^ -*-ln|r-r'| 
irK g 



■ 



where K g is the coupling constant of the Gaussian model (also called the stiffness constant). 
It varies continuously in the rough phase and takes the universal value ~ at roughening 
transitions. The interface free energies vanish in the rough phase as powerlaws. In a semi- 
infinite geometry M — > oo, rj m and rjm scale asymptotically as [f| 

mi = ~Y ^ 

where ( is the aspect ratio of the lattice constants in the spatial and time-like directions. 

We evaluate the interface free energies through the transfer matrix. Consider the transfer 
matrix for a square lattice rotated by 45° as shown in Fig. 4. In our units the aspect ratio 
is equal to ( = 2; one unit in time, a T is twice as big as the spatial unit a x . A height 
configuration (h , hi, ... , h N ) in a row is represented by a state vector \h , hi, ... , h N ). It 
is convenient to replace the height variables by step variables s« = hi — hi-i, with i = 
1, . . . , N. They take only the values 0, and ±1 due to the restricted solid-on-solid condition. 
The surface configuration in each row is therefore represented by \ho,s) where s stands for 
(si, S2, ■ ■ ■ , sn). The elements of the transfer matrix T are the Boltzmann weights associated 
with height configurations \ho,s) and \h' Q ,t) in successive rows. T is sparse, and can be 
expressed in terms of a product over local vertex-type scattering matrices, acting on the 
\ho,s) and \h' ,t) in successive rows and intermediate internal step variables u, defined in 
Fig. 4. 
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In the case of PBC's the step variables satisfy the conditions s i+ N = Sj and J2i=i s i — 
in all rows. The partition function, Zpbc = Tr Tp BC , and the free energy, are obtained from 
the largest eigenvalue e~ EpBC of Tpbc i n the M — ► oo limit. 

The transfer matrices Thi and Th2 for the horizontal BC's HI and H2 are easily defined. 
Only the conditions the step variables must satisfy change: In the case of HI, the step 
variables are again periodic, s i+ N = s i5 but with J2f=i s i — 2 in all rows. In the case of H2, 
the step variables are anti-periodic, s i+ N = — Sj, with no restriction in the value of J2iLi s i- 
The partition function in each case is given by Z a = Tr T^f and the free energy, in the 
M — > oo limit, is again obtained from the largest eigenvalue e~ Ea of T a (a = HI and H2). 
So the interface free energies are given by 

Vm = E m — -Epbc (6) 
Vm = E112 — Epbc ■ (7) 

The transfer matrices for the two vertical BC's are more intricate. They involve the 
symmetry properties Eq. (^|) and Eq. @ of the transfer matrix with PBC's. The translation 
invariance in the surface heights, Eq. (13), implies that Tpbc commutes with the symmetry 
operator 

P|M) = \h + 2,s) . (8) 

Therefore, it is useful to distinguish between two classes of surface states, {|e, s)} and {\o, s)}, 
i.e., all states with h even and odd, respectively. From the parity type symmetry property 
Eq. (H) it follows that T PB c commutes also with the operator R, defined by 

R|e, s) = |o, — s) 
R|o, s) = \e, — s) 

where — s stands for (— s 1; — s 2 , ■ ■ ■ ,—sn)- The transfer matrices for the vertical BC's can 
be expressed in term of Tpbc, P, and R as Zyi = Tr [Tppc P] and Z V2 = Tr [T^ C R]. 

To evaluate Zyi one needs to keep track of the height in the first column modulo 4. This 
makes this boundary condition less useful than its horizontal counter part HI, where we do 
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not need to keep track of the absolute height of the surface, and therefore can drop the h 
label of the surface configurations altogether. So we discard Zyi in the following analysis. 
On the other hand, Zy 2 is very useful. It can be written as 

i i 

where e~ EpBC< -^ (e~ £ pBcW) is the ith largest eigenvalue of T PB c in the R — +1 (—1) sector. 
By the R — +1 (—1) sector, we mean the set of state vectors which are eigenvectors of R 
with the eigenvalue +1 (—1). Unlike horizontal boundary conditions, rjy 2 depends on the 
entire eigenvalue spectra. However, in the thermodynamic limit, it can be approximated, 
up to the leading order, as 

1 e -ME PBC _ g-M£p BC 

where -Epbc (-^pbc) * s the largest eigenvalue of T PB c m the R — +1 (—1) sector. So the 
scaling behavior of i]y 2 is determined from the quantity 

m = Ep BC - E PBC , (10) 

i.e., the mass gap between the two i?-sectors. From the fact that ??V2 is finite in the flat 
phase and vanishes in the DOF and rough phases, it follows that this mass gap should be 
finite in the DOF phase and vanish in the flat phase. We will use both V2 and H2 to locate 
the preroughening phase boundary. 



III. NUMERICAL RESULTS AND SUMMARY 

The largest eigenvalues in each sector of the transfer matrix are obtained using the con- 
ventional iteration method. An arbitrary initial vector projects onto the largest eigenvector 
by applying the transfer matrix repeatedly. E PBC , E m , and En 2 are easily found by this 
method. Ep BC is obtained by choosing the initial vector in the R — — 1 sector. The state 
vector is (2 x 3^) dimensional for an semi-infinite strip of width N. The maximum strip 
width we can handle is iV = 12. 
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First, we focus on particular paths through the parameter space to illustrate the existence 
of the rough, DOF, and the flat phases. The FSS amplitude 

Sm(N) = ^ (11) 

of the HI type interface must converge to K g in the rough phase (see Eq. (|])). Numerical 
data of Shi (AO along the line Ea = 0.1 are presented in Fig. 5 (a). It shows that r/Hi scales 
as -k with continuously- varying amplitudes at small E B (high temperatures). The scaling 
breaks down at large Eb- This means that tja or rjs becomes nonzero. The roughening 
transition should take place when Shi (AO reaches the universal Kosterlitz-Thouless value 
|. This value is marked in Fig. 5 (a) by a dashed line, and indeed it crosses the numerical 
curve in the crossover region. So the numerical data in Fig. 5 (a) support that there is the 
rough phase at high temperatures, separated from the DOF or flat phase through a KT 
roughening transition. 

We present also the FSS amplitudes of the interface free energy r/H2 and the mass gap 
m, defined by 

S H2 (A0 = ^ (12) 
x(N) = ^ , (13) 

along the line Eb = 3.0 in Figs. 5 (b) and (c). Both quantities show crossing points. They 
signal the crossover between two regions. One where tja, the S^-step free energy, vanishes (at 
small Eb ) and one where it is finite (at large Eb)- This confirms the existence of the DOF 
phase and the PR transition, since these free energy gaps must be finite in the flat phase 
but converge exponentially to zero in the DOF phase. For some reason, the convergence for 
V2 is dramatically better than for H2. 

At KT type roughening transitions, the stiffness constant takes the universal value | . So 
we obtain a sequence of estimates of the roughening transition line by applying the condition 
Shi(AT) = | for each N. These are the roughening data points shown in Fig. 2. A sequence 
of estimates for the PR transition line between the flat and DOF phase can be obtained 



from the crossing points of x(N) and x(N — 2), and also those of Sm(N) and Sn2(N — 2). 
In Fig. 2 we show only the crossing points of V2 for N = Q,8, 10 and 12. (Those of H2 are 
much less convergent, see Fig. 5 (b) and (c)). 

The scaling theory of PR transitions, tells us that the critical fluctuations are described 
by the Gaussian model, but with K g greater than the universal KT value | of the roughening 
transition [0]]. We investigate this scaling behavior by studying the FSS amplitudes of 77H2 
and m. In the Gaussian model, Sn2(N), does not vary continuously, instead it takes the 
universal value n/4 (see Eq. @). On the other hand, the V2 type mass gap should scale 



as flTJ 



vr 2 C 

m = . (14) 

2K g N V ; 

This is related to the fact that at A = 0, the R = —1 and R = +1 sectors of Tpbc are 
equivalent apart from a phase factor e l7TSl attached to all step variables s\ at the seam (the 
first column of the lattice). 

In Fig. 6 (a) we present the FSS amplitude of the mass gap m. The vertical axis is scaled 
as 7i 2 (/(2mN), such that it represents K g , see Eq. (|H|). K g starts-off close to the universal 
value I in the neighbourhood of the roughening transition, at Eb — 0.8, and increases with 
Eb- This is in accordance with the assertion that the PR transition is described by the 
Gaussian model with continuously varying K g greater than |. The FSS behaviour of Sh2 
along the PR line is shown in Fig. 6 (b). The dashed line denotes the universal value | 
of the Gaussian model. The data at large Eb remain far from the universal value, though 
approach it. Like before, the convergence of this quantity is poor (see also Fig. 5(b)). 

The PR transition line in Fig. 2 seems to penetrate into the rough phase. But this does 
not mean that there is another transition inside the rough phase. The rough phase is a 
critical phase where the mass gap scales as 0(1/N) on either side of the crossing points. 
The presence of crossing points of V2 inside the rough phase represents only a turn around 
in the corrections to scaling amplitudes for the amplitude. 

The crossing of the two sets of lines in Fig. 2, the estimates for the roughening and PR 
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lines, confirms the existence of a multicritical point (Ea c , Eb c ) where the PR and roughening 
transition merge. A sequence of estimates (Ea c {N), Eb c (N)) for the multicritical point is 
obtained by solving the two conditions Shi (A) = | and x(N) = x(N — 2) simultaneously 
for each value of A = 6, 8, 10, and 12. These estimates are shown in Fig. 7. The arrows 
point towards power-law extrapolated values: 

E Ac = 0.41 ± 0.03 (15) 

Eb c = 0.89 ± 0.01 (16) 

In summary, we have investigated the phase transitions in a model system for Si or 
Ge (001) type crystal surfaces with an uniaxial structure that switches direction at each 
mono-atomic step. We obtained the phase diagram from a numerical FSS study of the 
transfer matrix spectra. It consists of flat, rough, and DOF phases. The unusual topological 
properties of the surface stabilize the DOF phase in the absence of step-step interactions, 
which are crucial for the stabilization of the DOF phase in conventional surfaces. The 
location of the multicritical point where the PR transition line merges to the roughening 
transition line is determined numerically, r c ~ 2.2. 

Specific crystals follow paths through Fig. 2 resembling straight lines as function of 
temperature, since the step energies are approximately constant. Our results shows that 
if the ratio r = Eb/Ea is greater than a critical value r c ~ 2.2, the flat unreconstructed 
crystal undergoes a PR transition into the DOF phase followed by a roughening transition. 
In Si (001) surfaces, the ratio between the step energies is larger, r ~ 3 We suggest 

therefore that Si (001) undergoes PR induced simultaneous deconstruction transition @. 
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FIGURES 

FIG. 1. Topology of Sa and Sb type step excitations on an unreconstructed Si (001) type 
surface. 

FIG. 2. The phase diagram of the model defined in Eq. (1). The roughening transition lines 
are obtained from Shi = § for N = 6 (n),8 (o), and 10 (A). Preroughening transition lines are 
obtained from the crossing points of x(N) and x(N — 2) for N = 6 (□), 8 (o), 10 (A), and 12 (y). 
The lines between the data points are guides to the eyes. 

FIG. 3. Step excitation type frustrations induced by the HI (a), H2 (b), VI (c), and V2 
(d) boundary conditions. The uniaxial direction in each domain of flat region is shown to help 
identifying the steps. 

FIG. 4. The transfer matrix set-up. The rows of a square lattice are rotated by 45°. The 
aspect ratio ( between the lattice constants in the horizontal and vertical directions is equal to 2. 
Height and step variables are defined on faces and bonds of the lattice, respectively. 

FIG. 5. Typical data of the FSS amplitude of the interface free energies r/ni (a), ?7H2 (b) and 
the mass gap m (c). Different symbols (□ for N = 6, o for N = 8, A for N = 10, and y for 
N = 12) are used to distinguish the strip widths N. The lines are guides to the eyes. 

FIG. 6. FSS amplitudes of m (a) and r)m (b) along the preroughening transition line. Different 
symbols (□ for N = 6, o for N = 8, A for N = 10, and y for N = 12) are used to distinguish the 
strip widths N. 

FIG. 7. Estimates (E Ac (N), E Bc (N)) for the location of the multicritical point for N = 6,8, 10, 
and 12. The extrapolated values are marked by arrows. 
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Figure 1 of Noh and den Nijs 
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Figure 2 of Noh and den Nijs 




Figure 3 of Noh and den Nijs 




Figure 4 of Noh and den Nijs 




Figure 5 of Noh and den Nijs 




Figure 6 of Noh and den Nijs 
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Figure 7 of Noh and den Nijs 



